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Abstract. This article presents several challenges to nuclear many-body theory and 
our understanding of the stability of nuclear matter. In order to achieve this, we 
present five different cases, starting with an idealized toy model. These cases expose 
problems that need to be understood in order to match recent advances in nuclear 
theory with current experimental programs in low-energy nuclear physics. 

In particular, we focus on our current understanding, or lack thereof, of many- 
body forces, and how they evolve as functions of the number of particles. We provide 
examples of discrepancies between theory and experiment and outline some selected 
perspectives for future research directions. 
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1. Introduction 

A central issue in basic nuclear physics research is to understand the hmits of stability 
of matter starting from its basic building blocks, either represented by effective degrees 
of freedom such as various hadrons or employing the underlying theory of the strong 
interaction, namely quantum chromo dynamics. To achieve this implies the development 
of a comprehensive description of all nuclei and their reactions, based on a strong 
interplay between theory and experiment. This interplay should match the research 
conducted at present and planned experimental facilities, where one of the aims is to 
study unstable and rare isotopes. These nuclei can convey crucial information about 
the stability of nuclear matter, but are difficult to produce experimentally since they 
can have extremely short lifetimes. They exhibit also unsual neutron-to-proton ratios 
that are very different from their stable counterparts. Furthermore, these rare nuclei lie 
at the heart of nucleosynthesis processes in the universe and are therefore an important 
component in the puzzle of matter generation in the universe. 

We do expect that these facilities will offer unprecendeted data on weakly bound 
systems and the limits of stability. To interpret such a wealth of experimental data and 
point to new experiments that can shed light on various properties of matter requires a 
reliable and predictive theory. If a theoretical model is capable of explaining a wealth of 
experimental data, one can thereafter analyze the results in terms of specific components 
of, say, the nuclear forces and extract simple physics pictures from complicated many- 
body systems. 

To better understand the rationale for this article, it is important to keep in 
mind the particularity of nuclear physics. Basic nuclear physics research, as conducted 
today, is very diverse in nature, with experimental facilities which include accelerators, 
reactors, and underground laboratories. This diversity refiects simply the rather complex 
nature of the nuclear forces among protons and neutrons. These generate a broad 
range and variety in the nuclear phenomena that we observe, from energy scales of 
several gigaelectronvolts (GeV) to a few kiloeletronvolts (keV). Nuclear physics is thus 
a classic example of what we would call multiscale physics. The many scales pose 
therefore a severe challenge to nuclear many-body theory and what are dubbed ab 
initio descriptions of nuclear systems. Examples of key physics issues which need to 
be addressed by nuclear theory are: 

• How do we derive the in-medium nucleon-nucleon interaction from basic principles? 

• How does the nuclear force depend on the proton-to-neutron ratio? 

• What are the limits for the existence of nuclei? 

• How can collective phenomena be explained from individual motion? 

J With ab initio we do mean methods which ahow us to solve exactly or within controlled 
approximations, either the non-relativistic Schrodinger's equation or the relativistic Dirac equation 
for many interacting particles. The input to these methods is a given Hamiltonian and relevant degrees 
of freedom such as neutrons and protons and various mesons. 
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• Can we understand shape transitions in nuclei? 

In order to deal with the above-mentioned problems from an ab initio standpoint, 
it is our firm belief that nuclear many-body theory needs to meet some specific criteria 
in order to be credible. Back in 2004 two of the present authors co-authored a preface 
to a collection of articles from a workshop on many-body theories, see Ref. [1] for more 
details. Seven specific requirements to theory were presented. These were: 

• A many-body theory should be fully microscopic and start with present two- and 
three-body interactions derived from, e.g., effective field theory; 

• The theory can be improved upon systematically, e.g., by inclusion of three-body 
interactions and more complicated correlations; 

• It allows for a description of both closed-shell systems and valence systems; 

• For nuclear systems where shell- model studies are the only feasible ones, viz., a 
small model space requiring an effective interaction, one should be able to derive 
effective two- and three-body equations and interactions for the nuclear shell model; 

• It is amenable to parallel computing; 

• It can be used to generate excited spectra for nuclei like where many shells are 
involved. (It is hard for the traditional shell model to go beyond one major shell. 
The inclusion of several shells may imply the need of complex effective interactions 
needed in studies of weakly bound systems); and 

• Finally, nuclear structure results should be used in marrying microscopic many- 
body results with reaction studies. This will be another hot topic of future ab 
initio research. 

Six years have elapsed since these requirements were presented. Recent advances in 
nuclear theory have made redundant most of the points listed above. They are actually 
included in many calculations, either fully or partly. On the other hand, in the last 
five years, we have witnessed considerable progress in nuclear theory. Of relevance 
to this article are the developments of effective field theories [2[ |^, lattice quantum 
chromodynamics calculations, with the possibility to constrain specific parameters of 
the nuclear forces [U El [6] and of many-body theories applied to nuclei, see for example 
Refs. [11 [HI [9], [To] . There are also several efforts to link standard ab initio methods with 
density functional theories, or more precisely, energy density functional theories, see 
Refs. [m 121 [13]. 

These developments have led us to formulate some key intellectual issues and further 
requirements we feel nuclear theory needs to address. The key issues are: 

(i) Is it possible to link lattice quantum chromodynamics calculations with effective 
field theories in order to better understand the nuclear forces to be used in a many- 
body theory? 

(ii) A nuclear force derived from effective field theory is normally constructed with a 
specific cutoff in energy or momentum space. Most interactions have a cutoff A 



Many-body interactions and nuclear structure 



4 



in the range A ~ 500-600 MeV. The question we pose here is how well do we 
understand the link between this cutoff and a specific Hilbert space (the so-called 
model space) used in a many-body calculation? 

(iii) The last question is crucially linked with our understanding of many-body forces 
and leads us to the next question. Do we understand how many-body forces evolve 
as we add more and more particles? Irrespective of whether our Hamiltonian 
contains say two- and three-body interactions, a truncated Hilbert space results 
in missing many-body correlations. To understand how these correlations evolve 
as a function of the number of particles is crucial in order to provide a predictive 
many-body theory. 

(iv) Finally, in order to deal with systems beyond closed-shell nuclei with or without 
some few valence nucleons, we need to link ab initio methods with density functional 
theory. 

With respect to requirements, we find it timely to request that a proper many-body 
theory should provide uncertainty quantifications. For most methods, this means to 
provide an estimate of the error due to the truncation made in the single-particle basis 
and the truncation made in limiting the number of possible excitations. The first point, 
as shown by Kvaal [11], can be rigorously proven for a chosen single-particle basis. The 
second point is more difficult and is normally justified a posteriori. As an example, 
in coupled-cluster calculations, a truncation is made in terms of various particle-hole 
excitations. Only specific sub-clusters of excitations are included to infinite order, see 
for example Ref . |15] . Whether this subset of excitations is sufficient or not can only be 
justified after the calculations have been performed. 

This article does not aim at answering the above questions. Rather, our goal is 
to raise the awareness about these issues since it is our belief that they can lead to a 
more predictive nuclear many-body theory. Our focus is on the third topic, and we 
illustrate the problems which can arise via a simple toy model in Section O Thereafter, 
we discuss four possible physics cases where the effect of missing many-body forces can 
be studied theoretically and benchmarked through existing and planned experiments. 
The physics cases are all linked with studies of isotopic chains of nuclei, with several 
closed-shell nuclei accessible to ab initio calculations. Our physics cases are discussed in 
Section 131 We start with the chain of oxygen isotopes, from ^^O to ^^O. The next cases 
deal with isotopes in the pf shell, namely the Ni isotopes (^^Ni, ^^Ni, ^^Ni, and ''^Ni) 
and the Ca isotopes ('^°Ca, ^^Ca, ^^Ca, ^^Ca, and ^"^Ca). We conclude with the chain of 
tin isotopes from ^"^^Sn to ^^°Sn. For isotopes like ^^O, ^^Ca, ^°Ca, and ^^°Sn, data on 
binding energies are missing. We try to give here a motivation why one should attempt 
to measure these nuclei. 

Our conclusions and perspectives are presented in the last section. 
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2. A simple toy model 

We start with a model described by a simple two-body Hamiltonian. This model catches 
the basic features concerning the connection between many-body forces and the size of 
the model space. In particular we will stress the link between a given effective Hilbert 
space and missing many-body effects. This section, although it forms the largest part 
of this article, conveys the problems which can arise in many-body calculations with 
truncated spaces and effective interactions. 

The model we present mimicks the effects seen in standard shell-model calculations 
with an effective interaction, either with a valence model space fL6\ or a so-called no-core 
model space [8]. 

2.1. Hamiltonian 

The Hamiltonian acting in the complete Hilbert space (usually infinite dimensional) 
consists of an unperturbed one-body part, Hq, and a perturbation V. The goal is 
to obtain an effective interaction Ves acting in a chosen model or valence space. By 
construction, it gives a set of eigenvalues which are identical to (a subset of) those of 
the complete problem. 

If we limit ourselves to, at most, two-body interactions, our Hamiltonian is then 
represented by the following operators 

H = ^{a\ho\(3)ala,3 + ^ ^{a(3\V\'y6)ala^^asa^, 

where aj^ and Oq,, etc. are standard fermion creation and annihilation operators, 
respectively, and a^jS represent all possible single-particle quantum numbers. The full 
single-particle space is defined by the completeness relation 1 = X^^i |q^)(q^|- In our 
calculations we will let the single-particle states |a) be eigenfunctions of the one-particle 
operator h^. 

The above Hamiltonian acts in turn on various many-body Slater determinants 
constructed from the single-basis defined by the one-body operator Hq. We can then 
diagonalize a two-body problem in a large space and project out via a similarity 
transformation an effective two-body interaction. The interaction acts in a much smaller 
set of single-particle states. The two-particle model space V is defined by an operator 

m 
a/3=l 

where we assume that m = dim{V) and the full space is defined by 

P + Q = i. 

Our specific model consists of N doubly degenerate and equally spaced single- 
particle levels labelled hj p = 1,2,... and spin a = ±1. These states are schematically 
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portrayed in Fig. [T] The first five single-particle levels define a possible model space 
indicated by the label V. 

We write the Hamiltonian as 

H = Ho + V, 

where 

^0 = ^^{p - l)al^ap„ 

pa 

and 

^ ^ XI ~^f^ (^al^al^ag_ar+ + h.c. j . 

pq pqr 

Here, Hq is the unperturbed Hamiltonian with a spacing between successive single- 
particle states given by ^, which we may set to a constant value ^ = 1 without loss 
of generality. The two-body operator V has two terms. The first term represents the 
pairing contribution and carries a constant strength g. (It easy to extend our model to 
include a state dependent interaction.) The indices cr = ± represent the two possible 
spin values. The first term of the interaction can only couple pairs and excites therefore 
only two particles at the time, as indicated by the rightmost four-particle state in Fig. [H 
There, a pair is excited to the state with p = 9. The second interaction term, carrying 
a constant strength /, acts between a set of particles with opposite spins and allows 
for the breaking of a pair or just to excite a single-particle state. The spin of a given 
single-particle state is not changed. This interaction can be interpreted as a particle- 
hole interaction if we label single-particle states within the model space as hole-states. 
The single-particle states outside the model space are then particle states. 

In our model we have kept both the interaction strength and the single-particle 
level as constants. In a realistic system like a nucleus, this is not the case; however, if a 
harmonic oscillator basis is used, as done in the no-core shell-model calculations [8], at 
least the single-particle basis mimicks the input to realistic calculations. 

2.2. Effective Hamiltonians 

We now consider a general A-body situation, of which our model is just an example, 
where a Hilbert space of finite dimension n is given along with an A-body Hamiltonian 
H with spectral decomposition given by 

n 
k=l 

where {\il)k)}^=i is an orthonormal basis of eigenvectors and where Ek are the 
corresponding eigenvalues. We choose the dimension n to be finite for simplicity, but 
the theory may be generalized to infinite dimensional settings where H has a purely 
discrete spectrum. 

The Hilbert space is divided into the model space V and its complement, denoted 
by Q. We assume again that m = dim('P). 
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Figure 1. Schematic plot of the possible single-particle levels with double degeneracy. 

The filled circles indicate occupied particle states while the empty circles represent 
vacant particle (hole) states. The spacing between each level p is constant in this 
picture. The first five single-particle levels define a possible model space, indicated by 
the label V. The remaining states span the excluded space Q. The first state to the 
left represents a possible ground state representation for a four-fermion system. In the 
second state to the left, one pair is broken by the interaction. All single- particle orbits 
belong to the model space. The two remaining four-particle states represents single- 
particle excitations to the excluded space, either by breaking two pairs or breaking one 
pair and exciting one pair (rightmost state). 
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The effective Hamiltonian ^^gfr is defined in the "P-space only, and by definition its 
eigenvalues are identical to the m eigenvalues of H. This is equivalent to Hcs being 
given by 

:= PnP 

= Pe-^He^P, 

where is assumed to obey the de-coupling equation 

QnP = 0. 

If the latter is satisfied, the P-space is easily seen to be invariant under "H, and since 
similarity transformations preserve eigenvalues, ^^efi is seen to have m eigenvalues of H. 

Without loss of generality, we assume that the eigenvalues Ek of H are arranged so 
that Hes, which is non-Hermitian in general, has the spectral decomposition 

m 

HeS = y^-Efc|0fc)(0fc|, 
k=l 

where {\(j)k)}^=i is a basis for the "P-space, and where {(t>k\(t>e) = Sk,e defines the bi- 
orthogonal basis {|0A;)}fcLi- 

The similarity transform operator exp(G) is, of course, not unique; E^, k = 1, ■ ■ ■ ,m 
can be chosen in many ways, and even if the effective eigenvector 10^) is chosen to be 
related to li'k), there is still great freedom of choice left. 

Assume that we have determined the eigenvalues E^, k = l,...,m that Hes 
should have. Two choices of the corresponding 10^) are common: The Bloch-Brandow 
choice, and the canonical Van Vleck choice, resulting in "the non-Hermitian" and "the 
Hermitian" effective Hamiltonians, respectively. For a discussion of these approaches 
see Refs. [I71[l8l[T9l[20l[2ll[221[23l[2l]. 

In the Bloch-Brandow scheme, the effective eigenvectors are simply chosen as 

\<pk) ■■=Pm, 

which gives meaning whenever P\ipk) defines a basis for P-space. In this case, G = u, 
where u = QuP, defined by 

Ljp\ijjk) := Q\ipk), fc = l,---,m. 

In contrast, the canonical Van Vleck effective Hamiltonian chooses a certain 
orthogonalization of {P\^pk)}T=l ^s effective eigenvectors |25j. In this case, G = 
arctanh(a) — u'^), which relates the two effective Hamiltonians to each other. The 
canonical effective interaction Hcs minimizes the quantity A defined by 

m 

Ai\xi),---,\Xm)):=Y.\\\Xk)-mr, (1) 

k=l 

where the minimum is taken with respect to all orthonormal sets of P-space vectors x- 
In fact, it can be taken as the definition [25]. The Bloch-Brandow effective eigenvectors, 
on the other hand, yield the global minimum of A. 



Many-body interactions and nuclear structure 9 

We have not yet specified which of tlie eigenvalues of H is to be reproduced by 
-f^eff- In general, we would like it to reproduce the ground state and the other lowest 
eigenstates of if if m > 1. We consider this question in some detail in the following 
sections. 

But first, we define the effective interaction \4ff as 

where [Hq, P] = is assumed. This is satisfied whenever the model space is spanned 
by Slater determinants being eigenvectors of Hq, such as in the model studied in this 
article. 

Finding Hes is equivalent to solving the original problem. In order to be useful, 
we need some sort of approximation scheme to find -^eff- It is common to use many- 
body perturbation techniques such as folded diagrams |T6], but as these suffer from 
convergence problems, the no-core shell model community [8] has developed the so-called 
sub-cluster approximation to V^e- This is a non-perturbative approach, and utilizes the 
many-body nature of H. The resulting effective interaction is often referred to as the 
Lee-Suzuki effective interaction [26] . 

In our model, V is defined by restricting the allowed levels accessible for the A 
particles under study. Thus, 

P(^) :={|(pi,ai)---(p^,a^)) : Pk < Np} , 

where Np < N is the number of levels accessible in the model space. This is the 
way model spaces in general are built up, simply restricting the single-particle orbitals 
accessible. It is by no means the only possible choice. On the other hand, this way of 
defining the model space has a very intuitive appeal, as it naturally leads to a view of 
^eff as a renormalization of H. It also gives a natural relation between model spaces 
for different A, which is absolutely necessary for the sub-cluster effective Hamiltonian 
to be meaningful. 

The effective Hamiltonian is seen to be an A-hodj operator in general, even though 
H itself may contain only two-body operators. Thus, V^ff can be written in its most 
general form as 

where ak = {pk,o'k) and u^Y'.. represent a specific matrix element. The approximation 
idea is then to obtain instead an a-body effective interaction V^ff^ where a < A, and 
view this as an approximation to V^^'' ■ This leads to 
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which is a much simpler operator, usually obtainable exactly by large-scale 
diagonalization of the a-body Hamiltonian. 

The remaining question is which eigenpairs of H^""^ should be reproduced by H^"^ , 
and which approximate eigenvectors should be used. There is no unique answer to this. 
The "best" answer would be for each problem to require a complete knowledge of the 
conserved observables of the many-body Hamiltonian [T^ [25]. 

On the other hand, if is a small perturbation, that is, we let V t— )■ AV^ and 
consider an adiabatic turning on by slowly increasing A, then it is natural to choose 
the eigenvalues developing adiabatically from A = 0. Indeed, V^^^ is then seen to be 
identical to a class of a-body terms in the perturbation series for the full vj^"^ to infinite 
order [231 1211 ■ The problem is, there is no way in general to decide which eigenvalues 
have developed adiabatically from A = 0, and we must resort to a heuristic procedure. 

Two alternatives present themselves as obvious candidates: Selecting the smallest 
eigenvalues, and selecting the eigenvalues whose eigenvectors have the largest overlap 
(ipklPlipk) with V. Both are equivalent for sufficiently small A, but the eigenvalues will 
cross in the presence of so-called intruder states for larger A [271 [281 [29] . Moreover, the 
presence of perhaps unknown constants of motion will make the selection by eigenvalue 
problematic, as exact crossings may lead us to select eigenpairs with P\ipk) = 0, which 
makes ill-defined. We therefore consider selection by model space overlap to be 
more robust in general. These considerations are, however, not important for our main 
conclusions. 

The non-Hermitian Bloch-Brandow effective Hamiltonian runs into problems in 
the sub-cluster approach since the interpretation of V^^'* as an interaction requires 
hermiticity when applied to an A-body problem - otherwise, eigenvalues will be complex 
unless a = A. 

See Ref. [25] for details on the algorithm for computing the effective interactions. 
2.3. Pure pairing interaction 

We first set / = 0, giving a pure pairing Hamiltonian. We study an A = 5 body 
problem with N = 8 levels and a model space consisting of the Np = 5 first levels. We 
compute vjg^ for a < 4 and compare their properties. We study here only the selection 
of eigenvectors using the overlap scheme. 

In Fig. [2] the error in the three first eigenvalues for each a are shown as a function of 
g G [—1,1]. For g > 0, a double-logarithmic plot reveals an almost perfect (yf^-behaviour 
of all errors. Larger values of a give smaller errors, as one would expect, but only by 
a constant factor. Thus, all the V^^'' seem to be equivalent to perturbation theory to 
second order in the strength g with respect to accuracy. This order is constant, even 
though the complexity of calculating V^g^ increases by orders of magnitude. Most of 
the physical correlations are thus well represented by a two-body effective interaction. 
This is expected since a pairing-type interaction favors strong two-particle clusters. The 
choice of a constant pairing strength enhances also these types of correlations. Three- 
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Figure 2. (Color online) Left: Error in the 3 lowest eigenvalues using V^g , a = 2 
(black), a = 3 (red), and a = 4 (blue) compared to exact eigenvalues (A = 5). Right: 
logarithmic plot for g > 0, showing almost perfect g'^-behavior. 



body and four-body clusters tend to be small. 

By investigating the remaining eigenvalues, we confirm that all effective eigenvalues 
behave in the described way, with only small variations. To assess the quality of the 
overlap selection scheme in this case, we graph the overlaps (ipklPl'ipk) of the selected 
and discarded eigenvectors for each a < 5 as a function of g in Fig. [31 Here, the selected 
overlaps are shown in black, while the discarded are colored red. We show only overlaps 
in the total spin S' = (a = 2,4) and 5 = 1/2 (a = 3, 5) channels for clarity. The 
selections for the P-space and the Q- space are well separated for all g, meaning that 
the selection scheme manages to follow the eigenpairs adiabatically from g = 0, where 
(ipklPl'ipk) is either zero or unity. Moreover, the view of the "P-space as being "effective" 
is sensible. However, for very large values of g, the graphs of the selected and discarded 
eigenvectors cross, meaning that the effective interaction view is broken. 

We conclude that the effective interaction V^^^ works very well in the pure pairing 
Hamiltonian long as the strength g is sufficiently small. However, there is not 

much to gain with respect to accuracy by going to a > 2 compared to a = 2. That is, 
a two-body effective interaction captures the relevant physics when only a pairing force 
is involved. In nuclear physics, the pairing interaction is rather strong, and the typical 
relation between the single-particle spacing ^ around the Fermi surface and the strength 
of the nuclear pairing interaction g is roughly — 5. In light nuclei like ^^O, 

the shell-gap is 11.25 MeV, and the average p-shell and sd-shell effective interactions 
are of the order of 1 — 5 MeV in absolute value [30]. In the region of the rare earth 
nuclei, typical single-particle spacings are of the order of some few hundred keV. The 
interaction matrix elements are of the same size in absolute value. Choosing thereby a 
parameter g slightly larger or smaller than one captures the essential physics produced 
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Figure 3. (Color online) Model space overlaps for a < 5, / = 0. Black lines are 
selected model space vectors' overlaps (V'fc|-P|V'fc)j while red lines are excluded space 
overlaps. Each line corresponds to one eigenvector of the a-body Hamiltonian. The 
■P-space selections and Q-space selections are clearly separated from each other. 



by a pairing interaction around the Fermi surface in nuclei. 
2.4- Pairing and particle-hole interaction 

We now turn on the particle- hole interaction, by setting f = ag with a > 0. In this case, 
the qualitative picture of the model space overlaps shown in Fig. [3] changes dramatically, 
even for very small values of a. In Fig. HJ we repeat the plots from Fig. [3] for a < 4 with 
a = 0.05. It is hard to see any systematic behaviour except for the a = 2 case. This 
goes to show that even tiny changes in an operator can give big changes in qualitative 
behaviour of its eigenvectors, even though the eigenvalues are only perturbed slightly. 
In the lower right plot of Fig. HI we have singled out the overlaps belonging to the 
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P-space overlaps, a=2 



P-space overlaps, a=3 





First 5 selected eigenvectors for a=3 




Figure 4. (Color online) Model space overlaps for a = 2,3, a = 0.05. Top figures: 
Black lines are selected model space overlaps, while red lines are excluded space 
overlaps. The a = picture (cf. top row of Fig. [3]) is radically changed. For a = 4 (not 
shown) the picture is even more complicated. These graphs show that it is difficult to 
choose a proper model space when the particle-hole interaction is turned on. Bottom 
figure: Singling out the five first selected overlaps for a = 3; the fifth (purple) shows 
a discontinuity. Note the vertical scale of the graphs. This shows that some of the 
chosen eigenpairs (typically the lowest) are chosen correctly. 



lowest eigenvalues for a = 3. It shows that there is still some ability left in the overlap 
selection scheme to select some eigenvalues adiabatically. The fifth overlap curve (and all 
the others not shown) has discontinuities, which would manifest itself as discontinuities 

This trend is general, as Table [1] shows. Here, the fraction d/m of the number of 
continuous overlaps d (counted visually) over model space dimension m for different 
a and a are shown for the subspace of lowest S = OorS* = 1/2 (for a = 3). 
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Table 1. Fraction d/m of "good" selections of eigenpairs for V^g , where m is the 
dimension of V. We limit the attention to states with quantum numbers S = = Sz = 
for a = 2,4, and S" = 1/2 = 5^ = 1/2 for a = 3. Also, / = 0.5g. The a = 2 case has 
a much larger fraction than the others, i.e., the overlap selection scheme works better. 
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Apparently, a > 2 gives a smaller fraction of continuous curves, meaning that following 
the eigenvalues adiabatically becomes increasingly difficult. The fact that the two-body 
case exhibits such good behavior is mainly due to the role played by the pairing force 
if the pairing force is stronger than the particle-hole contribution. With three or more 
active particles, our choice of particle-hole operator produces strong correlations between 
the model space and the excluded space. 

These considerations indicate that if any eigenvalues of H^'^ can be reliable, it will 
only be the lowest ones. Higher eigenvalues will almost certainly be without meaning. 
Of course, the matrix elements due to the "bad" selections (i.e., selections after a 
discontinuity in the curve) will also affect the lowest eigenvalues, but one can hope 
that these effects can be neglected. From Tabled] we thus expect V^g to perform better 
than values of a greater than two, but that only the lowest eigenvalues have a reasonable 
error behaviour. 

In Fig. [5] the eigenvalue error plots of Fig. |2] are repeated, but with a = 0.05. 
Clearly, a = 2 is able to produce sensible results in a wide range of g; the other V^^^ 
have errors that jump erratically. The picture gets worse when we try a = 0.5, shown 
in Fig. 12.41 In this case, not even the a = 2 case will give reliable results except for 
perhaps the ground state energy. We mention that in Fig. H] the a = 2 overlaps will 
cross at sufficiently large a, explaining thereby the breakdown for the a = 2 case as 
well. 

In summary, the particle-hole interaction changes the qualitative behaviour of the 
exact eigenvectors for the a-body problem drastically. It is clear that it becomes difficult 
to choose a model space V that fits well with the model space overlaps shown in Fig. |H 
such that it becomes possible to follow the eigenvalues adiabatically as a function of 
g. The two-body effective interaction, however, is surprisingly well behaved for a large 
range of interaction strenghts. 

We have seen that using using effective interactions with a > 2 does not seem to 
yield extra accuracy in general compared to the a = 2 case. First of all, such vj^^ 
interactions are much harder to compute. Secondly, we have observed only a constant 
factor gain in accuracy in the pairing case, in addition to obvious problems when the 
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Figure 5. (Color online) A few errors in eigenvalues using v}g\ f = 0.05g. Lower 
right: corresponding exact eigenvalues for A — 5. Each curve is colored corresponding 
to the exact eigenvalue curves. The a = 2 is well-behaved while the a = 3 case and the 
a = 4 case give unreliable eigenvalues. 



particle- hole interaction is present. 

However, with a strong particle-hole interaction, intruder-states can start playing 
a major role, and effective many-body interactions with a > 2 are not capable of 
reproducing the low-lying levels. The crossings of the curves in, for example. Fig. H] 
are due to intruder states, and these are totally absent in integrable systems in contrast 
to the apparent abundance seen here. Indeed, all two-particle systems described by 
central force interactions and external harmonic oscillator potentials are classically 
integrable. There are several examples of nuclear systems where intruder states play a 
major role, with perhaps the so-called island of inversion for nuclei with mass number 
A ~ 31 — 33 as one of the more popular mass regions studied recently (see for example 
the experimental results on the /3-decay of ^'^Mg in Ref. [31]). Both the parent and 
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(Color online) A few errors in eigenvalues using V^^"*, / = 0.5g. Lower right: 
corresponding exact eigenvalues for A = 5. Each curve is colored corresponding to 
the exact eigenvalue curves. The a = 2 is now less well behaved. An erratic behaviour 
is, in all cases, amplified compared to the / = 0.05^^ case. Compare also the error-axis 
range with Fig. [5l 



the daughter (^^Al) nuclei reveal intruder-state configurations of two-particle-two-hole 
character in the lowest excited states, in addition to proposed admixtures of one-particle- 
one-hole and three-particle-three-hole configurations for the ground state of ^^Mg. 

From our results, we may therefore infer that when intruder states are present, 
special attention must be paid to the construction of an effective interaction, as pointed 
out almost four decades ago by Schucan and Weidenmiiller [271 EH] • 

To link this discussion with the topic of missing many-body correlations, we can 
view the presence of intruder states as an example of a model space which is too small. 
This is in turn reflected in missing many-body correlations beyond those which can be 
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produced by say a given sub-cluster effective interaction. Reducing the model space 
thus produces missing many-body correlations, which, depending on the strength of the 
interaction, can produce large differences between the exact results and those obtained 
with a sub-cluster effective interaction. 

We can therefore summarize by saying that an effective sub-cluster Hamiltonian will 
always produce missing many-body physics. The size of these contributions can normally 
only be determined a posteriori, that is after a calculation has been performed. To be 
able to estimate the role of these effects is important for nuclear many-body theory, 
as it provides us with a sound error estimate. The hope is obviously that their effect 
is negligible, although the calculations of Ref. [32] indicate small but non-negligible 
four-body contributions for ^He when a three-body Hamiltonian that reproduces the 
experimental binding energy of '^H is employed. When one employs nuclear forces 
derived from effective field theory, many-body terms such as three-body interaction 
arise naturally [2]. There are now clear indications from several calculations (see for 
example Refs. [3], El El |32l |33l |3ll |35l [36]) that at least three-body interactions have to 
be included in nuclear many-body calculations. Whether four-body, or more complicated 
many-body, interactions are important or not for an accurate description of nuclear data 
is an open issue. From a practical and computational point of view, we obviously favour 
negligible four or higher many-body interactions. 

To understand how many-body interactions develop as one adds more and more 
particles is thus an open and unsettled problem in nuclear many-body physics (or many- 
body physics in general if the effective Hilbert space is too small). These interactions 
are not observables; however, their role can be extracted from calculations with a given 
effective sub-cluster interaction. 

To make this clear, assume that a given three-body Hamiltonian derived from 
effective field theory has been used to study the chain of oxygen isotopes. With this 
Hamiltonian, we perform then precise many-body calculations for all isotopes from ^^O 
to ^^O. If the errors of our calculations are negligible, we can infer that these are the 
results with this specific effective Hamiltonian. The discrepancy between theory and 
experiment can then be used to extract the role of missing many-body forces as a 
function of the number of nucleons. The missing many-body physics depends obviously 
on the employed effective Hamiltonian. 

Experiment or simple models provide us therefore with the benchmark our favourite 
theory has to reproduce. If the experimental data are reproduced to within given 
uncertainties, one can start analyzing the properties of nuclei in terms of various 
components of the nuclear forces. This can allow for the extraction of simple physical 
mechanisms, as done recently by Otsuka et al [37] . 

This coupling between experiment and precise many-body calculations is what we 
will sketch in the next section, starting with the oxygen isotopes. 
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3. Physics cases 

The nuclear forces derived from effective field theory are obviously much more 
complicated than the above simple model. First of all, our simple model contains only 
a two-body Hamiltonian. It allowed for a numerically exact diagonalization for a five- 
particle system in a Hilbert space consisting of eight doubly degenerate single-particle 
states. An effective sub-cluster interaction defined for a smaller space (typically the 
five lowest doubly degenerate single-particle states) gave results which were close to the 
exact ones if the interaction was weak and a relevant model space was employed. In 
particular, with a particle-hole interaction, one can easily get a strong admixture from 
states outside the chosen model space. In that case, the various sub-cluster effective 
interactions failed in reproducing the low-lying eigenvalues. 

The nuclear interactions derived from effective field theories [2] already contain 
at the so-called N^LO level, three-body interactions. At the N^LO level, four-body 
interactions appear. These interactions are normally constructed for a given cutoff 
A. In principle, if one uses an interaction at the N^LO level, one should also include 
four-body interactions. Normally, only three-body interactions at the N^LO level are 
introduced, together with a two-body interaction that fits two-particle data (scattering 
data and bound state properties) at the N^LO level. 

This means that when employed in a many-body context with more than three 
particles, our Hamiltonian lacks some many-body correlations. Hopefully these are 
small. A comparison with data, if our theory produces converged results at a given level 
of many-body physics, may then reveal their importance. 

In this sense, a given nuclear interaction computed with effective field theory at 
a given N"LO level, should contain all possible many-body interactions. At the N^LO 
level, one should include one-body, two-body, three-body, and four-body interactions. 
Omitting some of these leads to a sub-cluster effective Hamiltonian. In a certain sense, 
this parallels the discussion from the previous section. However, in the nuclear physics 
case, it is important to keep in mind that these interactions have been derived with a 
specific energy cutoff. The energy cutoff defines our effective Hilbert space, or model 
space. If the model space defined for an actual many-body calculation does not include 
all possible excitations within this model space, one can easily end up with the type of 
missing many-body correlations discussed in the previous section. 

To make our scheme more explicit, assume that we are planning a calculation of the 
oxygen isotopes with an effective field theory interaction that employs a cutoff A = 600 
MeV. Assume also that our favourite single-particle basis is the harmonic oscillator, 
with single-particle energies Eni = huj{2n + I + 3/2), with u the oscillator frequency, 
n = 0, 1, 2, . . . being the number of nodes and / the single-particle orbital momentum. 
The oscillator length b is defined as 



We define p = 0,l,2,...,P with P = 2n + / as the quantum number p of the highest- 




Many-body interactions and nuclear structure 



19 



filled level. The level labelled p can accommodate (j9 + + 2) fermions, with a 
spin degeneracy of two for every single-particle state taken into account. For a given 
maximum value of P = 2nmax + ^max, we have a total of 



single-particle states. 

The cutoff A defines the maximum excitation energy a system of A nucleons can 
have. The largest single-particle excitation energy (corresponding to possible one- 
particle-one-hole correlations) is then 



The value of huj can be extracted from the mean-squared radius of a given nucleus. One 
can show that this results in [38] 



with ro ~ 1 fm. Setting A = 16 and A = 600 MeV results in P 42. The largest 
possible value for n is then nmax ~ 21, or 22 major shells. With P = 42, the total 
number of single-particle states in this model space is 28380! 

For ^^O, this means that we have to distribute eight protons and eight neutrons in 
28380 single-particle states, respectively. The total number of Slater determinants, with 
no restrictions on energy excitations, is 



Any direct diagonalization method in such a huge basis is simply impossible. One 
possible approach is to introduce a smaller model space with a pertinent effective 
interaction, similar to what was discussed in the previous section. Depending on the 
size of the model space and the strength of the interaction, this can lead to uncontrolled 
many-body correlations, with effects similar to what was discussed within the framework 
of the simple model from Section O This has been the philosophy of the no-core shell- 
model approach [8j or many-body perturbation theory [16]. Another alternative is to 
employ methods which allow for systematic inclusions to all orders of specific sub-clusters 
of correlations within the full model space defined by A. This reduces the computational 
effort considerably and accounts, at the same time, for most of the relevant degrees of 
freedom. The Coupled cluster method [IHl [391 SO] and Green's function theory [lO] allow 
for such systematic expansions starting from a given Hamiltonian. 

Coupled-cluster theory has been particularly successful in both quantum chemistry 
and nuclear physics. It fulfills basically all the requirements we listed in Section [TJ 
Furthermore, with a given truncated single-particle basis (our effective Hilbert space) we 
can extract precise error estimates on the energy [lH [H] . It is a topic of current research 




p=0 



A = hwP = hw{2nms.x + /max)- 
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to be able to predict and understand the error due to truncations in the number of many- 
body correlations. As an example, most coupled-cluster calculations seldom go beyond 
correlations of the three-particle-three-hole type, the so-called triples correlations. An 
estimate of the error made due to the omission of four-particle-four-hole correlations 
(and higher-order correlations) would therefore be very useful. That would lead to a 
much more predictive theory. 

3.1. Brief review of coupled cluster theory 

Our many-body method of choice is thus coupled-cluster theory, see Refs. [T5l|39lll0] for 
more details. The coupled-cluster method fulfills Goldstone's linked-cluster theorem and 
therefore yields size-extensive results, i.e., the error due to the truncation is linear in the 
mass number A of the nucleus under consideration. Size extensivity is an important issue 
when approximate solutions to all but the lightest nuclei are sought [I5l |12] . Second, 
the computational effort scales gently (i.e., polynomial) with increasing dimension of 
the single-particle basis and the mass number A. The method has met benchmarks in 
light nuclei [33lll3]. 

Our algorithm is as follows: 

(i) We start with a nuclear force from effective field theory at a given N"LO order. 
In the calculations we present here, we include only a nucleon-nucleon interaction 
to order N'^LO. We neglect three-nucleon forces since their application within the 
coupled-cluster method is still hmited to smaller model spaces [M]. This gives rise 
to one source of systematic error and most likely the largest uncertainty in our 
results. So-called power-counting estimates from chiral effective field theory result 
in an uncertainty of about 2 MeV per nucleon, see Ref. [9] for further details. 

(ii) We compute ground state properties and excited states by including one-particle- 
one-hole, two-particle-two-hole, and selected three-particle-three-hole correlations. 
These correlations are called singles, doubles, and triples, respectively. We 
consider corrections due to triples excitations within the so-called ACCSD(T) 
approximation, see for example [9] for more details. Comparison of a calculation 
which includes only one-particle-one-hole and two-particle-two-hole correlations 
(CCSD), see Ref. [9], shows that triples corrections account for 10-15%. Similar 
ratios are found in coupled-cluster calculations of atoms and molecules, and 
experience in quantum chemistry (see for example Ref. [15]) suggests that the 
truncation of the cluster amplitudes beyond the triples corrections introduces an 
error of a few percent. 

(iii) The actual single-particle space used is defined close to the effective field theory 
model space. Our error here is smaller than the error made in the previous two 
points. The final results are more or less converged as a function of the size of the 
model space. A precise error estimate for a given single-particle basis was presented 
within a configuration interaction approach by Kvaal [H]. Similar error estimates 
apply to coupled-cluster theory as well [H]. 
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Within the choice of model space, most of our expectation values are practically 
converged. This is, for example, reflected in the fact that the ground state energy is 
almost independent of the chosen harmonic oscillator energy (for demonstrations, see 
for example Refs. [9l |33l SO]). From this, we can infer that our specific many-body 
approximation yields probably the optimal result which can be achieved at this specific 
level of many-body physics. Eventual discrepancies with experiment can then most 
likely be attributed to missing many-body correlations. A typical case would be the 
lack of three-body interactions from effective field theory. 

3.2. Oxygen isotopes 

The oxygen isotopes are the heaviest isotopes for which the drip line is well established. 
There are large experimental campaigns worldwide [HI HSl |16] which aim at uncovering 
the properties of the oxygen isotopes, both at or close to the drip line and beyond. 
Two out of four stable even-even isotopes exhibit a doubly magic nature, namely ^^O 
{Z=8,N=U) HZ] and ^^O {Z=8, N=16) HH |19]. The structure of these two doubly 
magic nuclei is assumed to be governed by the evolution of the lsi/2 and 0^5/2 one- 
quasiparticle states. The isotopes ^s-zsq ^^-^.^ believed to be unstable towards neutron 
emission, even though ^^O is a doubly magic nucleus within the standard shell-model 
picture. 

Of interest to us is the fact that we can perform ab initio coupled-cluster calculations 
for all assumed closed-shell nuclei of this isotopic chain, that is, oxygen isotopes with 
A = 16, A = 22, A = 24, and A = 28. Furthermore, we can also compute the Ail 
nuclei such as ^^F, ^^O, and ^^N. This means further that we can extract the isospin 
dependence of say spin-orbit partners. In Ref. [9] we considered the nuclei 16,22, 24,28q ^^^^ 
computed their ground-state energies within the ACCSD(T) approximation for chiral 
interactions with cutoffs of = 500 MeV and A^ = 600 MeV, respectively. Figure 
[6] shows the ground state energies relative to ^20. At a cutoff A^ = 500 MeV, we find 
that ^^O is bound by about 2.7 MeV with respect to ^'^O. However, the situation is 
reversed at the higher cutoff A^ = 600 MeV, and the difference is about -1.3 MeV. With 
the uncertainties due to missing three-body interactions, it is presently not possible 
to reach a conclusion regarding the existence of ^^O. However, using interactions from 
chiral effective field theory, we see that the stability of ^^O depends mainly on the 
contributions of the three-nucleon force (and probably more complicated many-body 
forces), and that even small contributions can tip the balance in either direction. A 
mass measurement of ^^O is clearly needed. 

The above results show that we do not yet understand fully how many-body forces 
evolve as we add more nucleons. The hope is that the inclusion of a three-body 
interaction from effective field theory reduces the dependence on the cutoff at the two- 
body level. This was nicely demonstrated by Jurgenson et al [32] in a recent study of 
^H and ^He. For '^He, there were indications that a small four-body interaction may be 
needed, as also demonstrated by Epelbaum et a/ [3]. Our ab initio calculations indicate 
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Figure 6. (Color online) Ground-state energies of neutron-rich oxygen isotopes "^O 
relative to ^^O for chiral interactions with two different cutoffs A^. Taken from Ref. [9]. 



also that the recent resuhs from phenomenological shell-model approaches regarding 
the unbound character of ^^O should be viewed with caution. The combination of 
three-nucleon forces, the proximity of the continuum, and the isospin dependence are a 
challenge for reliable theoretical predictions, see also Refs. |50l l5T] . 

Summarizing, we can claim that there is a posteriori evidence of the need of a three- 
body interaction or more complicated interactions. We hope that with the inclusion of 
three-body forces, the difference seen for different chiral interactions should become 
small. 

Furthermore, if the differences between theory and experiment are small, theory 
can be used to extract simple information on say various components of the nuclear 
force. We can, for example, study the evolution of many-body forces as we increase 
the number of valence nucleons. One possibility is to extract the A-dependence and 
the isospin dependence of such correlations as we move towards the drip line. Another 
interesting quantity which can be extracted is the isospin and A-dependence of the 
spin-orbit force. This will allow us to answer and understand which mechanisms change 
the single-particle fields close to the drip line. For these mass regions and the oxygen 
isotopes, experimental data on ^^O are very important. Similarly, for spin-orbit partners, 
measurements of closed-shell nuclei with one particle added or removed beyond A = 24 
will allow us to benchmark theory with experiment and perhaps understand which parts 
of the nuclear force are important close to the drip line. 



3.3. Calcium and nickel isotopes 

Doubly magic nuclei are particularly important and closed-shell nuclei like ^^Ni, 
^°°Sn and ^^^Sn have been the focus of several experiments during the last several 
years [521 ESI EH |55l [561 [HH [HHl [591 [60] • Their structure provides important information 
on theoretical interpretations and our basic understanding of matter. In particular. 
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recent experiments [531 EH ESI ESI EZ] have aimed at extracting information about 
single-particle degrees of freedom in the vicinity of ^^Ni. Experimental information from 
single-nucleon transfer reactions and magnetic moments [521 |53l EH ESI EZ| can be used 
to extract and interpret complicated many-body wave functions in terms of effective 
single-particle degrees of freedom. Transfer reactions provide, for example, information 
about the angular distributions, the excitation energies, and the spectroscopic factors 
of possible single-particles states. 

Much of the philosophy exposed in the previous subsection can be repeated for 
studies of calcium and nickel isotopes. The interesting feature here is that we probe 
systems with more nucleons. The chain of nickel isotopes exhibits four possible closed- 
shell nuclei, namely ^^Ni, ^^Ni, ^*Ni, and ''^Ni, with a difference of 30 neutrons from ^^Ni 
to '''^Ni. Similarly, for calcium isotopes in the /p-shell, we have 20 neutrons that can 
be added to ^°Ca and five possible closed-shell nuclei, ^°Ca, "^^Ca, ^^Ca, ^^Ca and ^°Ca. 
The binding energies of ^^Ca and ^°Ca are presently not available. Moreover, compared 
with the oxygen isotopes with at most 12 nucleons in the sci-shell outside ^^O, different 
single-particle degrees of freedom are probed. Comparing theory and experiment can 
again provide important information about spin-orbit partners close to the drip line, 
with information on their density and isospin dependence as well. 

As stated in the previous subsection, we are in the position where such nuclei can be 
calculated rather accurately with a given two-body Hamiltonian. This is demonstrated 
in Refs. [TOl HQ]. We are not limited to closed-shell systems but can also compute 
ground states and excited states of A ± 1 {A is the closed-shell nucleus) nucleons 
rather precisely. It means that a nucleus like ™Cu can be accessed with ab initio 
methods like coupled cluster theory. Single-particle properties like magnetic moments 
and spectroscopic factors can then give a measure of how good a closed nucleus ^^Ni is. 
Experimental data which probe single-particle properties on A ± 1 systems close to the 
dripline, will therefore provide important benchmarks for theory. 



3.4- Ab initio studies of nuclei around mass A = 100 

We end this section by stating that the analysis which can be performed for /p-shell 
and sc?-shell nuclei can be extended to the region of the tin isotopes as well, with both 
^"^^Sn and ^^^Sn as two very important nuclei for our understanding of the stability of 
matter, see for example the recent works of Refs. [581 EHl EOl [61]. In this case we expect 
to be able to run similar calculations within the next two to three years for nuclei like 
^"^^Sn, ^^^Sn, and ^^°Sn together with the corresponding A±l nuclei. We can then test 
the development of many-body forces for an even larger chain of isotopes and provide 
theoretical benchmarks for nuclei near ^^°Sn, of great importance for the understanding 
of the r-process, a nucleosynthetic process responsible for the production of around half 
of the heavy elements. A recent experiment on ^^^Sn shows clear evidence that this 
nucleus is an extremely good closed-shell nucleus [61]. This has, in turn, consequences 
for our understanding of quasi-particle properties in this mass region. We mention also 
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here that ^^''Sn is the last reported neutron-rich isotope (with half-hfe). Our final aim is 
to provide reliable predictions for all possibly closed-shell nuclei, from A = 4to A = 208. 

4. Conclusions 

The aim of this article has been to shed light on our understanding of many-body 
correlations in nuclei. Since all theoretical calculations involve effective Hamiltonians 
and effective Hilbert spaces, it is crucial to have a handle on the role many-body 
correlations play in a many-body system. This means that a sound theory should provide 
error estimates on the importance of neglected many-body effects. To understand these 
and develop sound error estimates is mandatory if one wants to have a predictive theory. 

We are now in a position where fairly precise results can be obtained for several 
closed-shell nuclei with a given two-body Hamiltonian. However, since most two-body 
interactions used nowadays are based on chiral effective field theory, it means that 
three-body and, more complicated, many-body interactions arise at different chiral 
orders [2]. The neglect of such higher-order terms leads to different predictions for 
two-body Hamiltonians. Experiment can then be used to understand how important 
these neglected degrees of freedom are. Stated in a more philosophical way, we seem to 
be doing pretty well at 'doing the problem right' (verification); but we still struggle with 
'validation' (doing the right problem). Having now developed the tool set to rigorously 
solve the nuclear many-body problem accurately, we can use these tools to more fully 
investigate the nuclear interaction. 

The hope is that the addition of three-body interactions will produce results which 
are more or less independent of the cutoff used in chiral effective field theory. This means, 
in turn that if our theoretical calculations with a three-body Hamiltonian reproduce, for 
example, experimental binding energies to high accuracy (theoretical and compared with 
experiment), we can then start analyzing properties like single-particle energies close to 
the drip line or the A-dependence of missing many-body correlations. Our formalism 
allows us now to compute closed-shell nuclei along an isotopic chain and extract single- 
particle energies and properties. For the chains of oxygen, calcium, nickel, and tin 
isotopes, we can therefore provide important information about nuclear structure at 
the drip lines of these chains. But in order to do this, one needs experimental data 
on properties like the binding energies of ^^O or ^^Ca. We feel that the enormous 
progress which has taken place in the last few years in nuclear theory can really lead to 
a predictive and reliable approach to nuclear many-body systems. In this endeavor, the 
close ties between theory and experiment are crucial in order to understand properly 
the limits of stability of matter. 
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